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ABSTRACT 

In a recent paper [ Phys . Fluids A2: 1678- 1684, 1990], the authors compared the perfor- 
mance of a variety of turbulence models including the K — e model and the second-order 
closure model derived by Yakhot and Orszag based on Renormalization Group (RNG) meth- 
ods. The performance of these RNG models in homogeneous turbulent shear flow was found 
to be quite poor, apparently due to the value of the constant C e \ in the modeled dissipation 
rate equation which was substantially lower than its traditional value. However, recently a 
correction has been made in the RNG based calculation of G c \- It is shown herein that with 
the new value of C e i, the performance of the RNG K — e model is substantially improved. 
On the other hand, while the predictions of the revised RNG second-order closure model are 
better, some lingering problems still remain which can be easily remedied by the addition of 
higher-order terms. 


*Thia research was supported by the National Aeronautics and Space Administration under NASA Con- 
tract No. NAS 1-18605 while the first author was in residence at the Institute for Computer Applications in 
Science and Engineering (ICASE), NASA Langley Research Center, Hampton, VA 23665. 
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A comparative study of the performance of nine independent turbulence models in ro- 
tating homogeneous shear flow was recently reported by Speziale et al. 1 Two of the models 
considered consisted of the K — c model and second-order closure model derived by Yakhot 
and Orszag 2 using Renormalization Group (RNG) methods. It was rather surprising how 
poorly the RNG models performed in homogeneous shear flow relative to the older, em- 
pirically based models of the same general type. The origin of the deficient predictions of 
the RNG models appeared to be largely due to the rather low value of the constant C t \ in 
the modeled dissipation rate equation; the RNG value of C t \ was 1.063 in contrast to the 
more traditional value of C e i = 1.44. However, a recent re-examination of the RNG based 
calculation of C e \ by Yakhot and Smith 3 has led to a correction - the new value of C e 1 is 
1.42. Some minor changes in the values of other constants in the RNG K — e model were also 
made. 3 In light of these changes, it would be desirable to set the record straight in regard 
to what these Renormalization Group models now predict for homogeneous shear flow - a 
critical test case used to evaluate the performance of models. This establishes the motivation 
for the present paper. 

In the RNG K - e model, the Reynolds stress tensor = t i[u'j (given that u- is the 
fluctuating velocity and an overbar represents an ensemble mean) is modeled as follows: 2,3 



where K = is the turbulent kinetic energy, e = vdu'Jdxj du'Jdxj is the turbulent dis- 

sipation rate, U{ is the mean velocity, and is a dimensionless constant which is calculated 
to be 0.085. In homogeneous turbulence, the turbulent kinetic energy is a solution of the 
transport equation 

(2) 

which is exact. The turbulent dissipation rate is obtained from the RNG derived transport 
equation 


i=-C. i-Ti 


dui 


' tX K Tij dxj W2 'tf 

where C e \ = 1.42 and C e i = 1.68 according to the recent calculations of Yakhot and Smith. 3 
These new values constitute a correction to the earlier values of C e i — 1.063 and C , j = 1.72 
reported by Yakhot and Orszag. 2 An additional production term was also uncovered by 


( 3 ) 


1 


Yakhot and Smith 3 which they were unable to close. However, an order of magnitude 
analysis 3 indicated that this term is small unless there are large strain rates - a case which 
will not be considered herein. Hence, we will neglect this additional term in the present 
study. For the RNG second-order closure model, the eddy viscosity model (1) is replaced 
with a Reynolds stress transport model of the form 1 



where C\ and C 2 are constants that are calculated to be 1.59 and 2/15, respectively. Some 
clarifications are needed concerning the origin of this model which has not been published 
and was obtained from a private communication with V. Yakhot. We have come to learn 
that this was not intended to be a final model, but rather was the result of a low-order 
calculation of the pressure- strain correlation whose purpose was to merely demonstrate that 
the Rotta term - with a coefficient C\ close to the well accepted value of 1.5 - could be 
formally obtained from RNG. Hence, the results predicted by this preliminary model should 
be judged accordingly. 

In homogeneous shear flow, an initially isotropic turbulence where 

2 

T\j = 3^0^ 1 £ = £0 ( 5 ) 


at time / = 0 is subjected to a constant shear rate S with the corresponding mean velocity 
gradient tensor 


dui 

dxj 


0 5 0 \ 
0 0 0 
0 0 0 / 


( 6 ) 


In Figure 1, the time evolution of the turbulent kinetic energy (where K* = K/Kq and 
t* = St) predicted by the new RNG K — e model is compared with the large-eddy simulation 
of Bardina et a/. 4 for an initial condition of Eq/SKq = 0.296. The predictions of the old 
version of the RNG K — e model (where C M = 0.0837, G c \ = 1.063, and C e 2 = 1.72) as well 
as the standard K — e model (where = 0.09, C e \ — 1.44 and G e 2 = 1.92) are also shown 
in Figure 1. It is clear from these results that the revised RNG K — e model does the best 
overall job in reproducing the growth rate of the numerical experiment on homogeneous shear 
flow. Analytically, it can be shown why this is the case. From a straightforward calculation, 
it can be shown that the turbulent kinetic energy and dissipation rate grow exponentially in 


2 



homogeneous shear flow as follows: 1,5 

K* ~ exp(A t*), e* ~ exp(Ai*) 
where the dimensionless growth rate A is given by 


A = 


- c A y 1 

(C tl - l)(C t2 - 1). 


1/2 


( 7 ) 


Hence, the growth rate becomes singular when C e \ = 1 — a state of affairs that explains 
why the old version of the RNG K - e model, with C e i = 1-063, overpredicted the growth 
rate of the turbulent kinetic energy by such a wide margin. The new version of the RNG 
K — e model predicts a growth rate of 


A = 0.142 


which is extremely close to the range of values obtained from physical and numerical 
experiments. 6,7 On the other hand, the standard K — c model predicts the somewhat high 
value of A = 0.226 which explains why this model overpredicts the LES data for K* as shown 
in Figure 1. A more complete set of the equilibrium values predicted by these different ver- 
sions of the K — £ model will be provided later. 

In Figure 2, the time evolution of the turbulent kinetic energy predicted by the revised 
RNG second-order closure model is compared with the large-eddy simulation of Bardina et 
al. 4 as well as with the predictions of the earlier version of the model and the Launder, Reece, 
and Rodi 8 (LRR) model. The new version of the RNG model does yield better predictions 
than the older version of the model since the previous value of C e \ = 1.063 was too close 
to C e i = 1 which constitutes a bifurcation point of the dissipation rate transport equation 
as shown by Speziale. 9 However, there are still problems with the model which gives rise to 
points of inflection in the time evolution of K* - a feature that makes it inferior to other 
second-order closure models such as the Launder, Reece, and Rodi model. The origin of 
this problem appears to be tied to the modeling of the pressure strain correlation. In the 
Launder, Reece, and Rodi model, the pressure strain correlation = p^du'Jdxj + du'Jdxi) 


is modeled as follows 


n v 


— 2Ciebij -f- 2C 2 K^ tJ + C$K 



+ bjk^ik — o 



+C<K{b ik W jk + b jk W ik ) 


( 8 ) 
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where 


( 9 ) 


* J 2 \c?Xj dxi ) ’ * J 2 \(9xj 9xj / 

Ji ' = ( T « “ i**«) / 2ir (i°) 

(in the simplified form of the Launder, Reece, and Rodi model, C\ = 1.8, C% = 0.4 and 
( 7 3 = C4 = 1.2). This model satisfies two important consistency conditions: (a) the constant 
C2 is equal to 2/5 - a result that follows from simple symmetry conditions for ILj as well 
as from Rapid Distortion Theory, 8 and (b) it represents a formal expansion of IX,y to 0 (b) 
in the anisotropy tensor. On the other hand, for this preliminary RNG second-order closure 
model we have 

Tlij = —2C\ebij + 2 C 2 K^ij ( 11 ) 

where C\ = 1.59 and C 2 = jg. This model is not complete to 0 (b) in the rapid pressure- 
strain term and violates the important symmetry constraint of C 2 — 2 / 5 . The fixed points 
that the resulting nonlinear ODE’s for these second-order closure models give rise to in 
homogeneous shear flow are of the focus type. 5 Significant deviations of C 2 from a value 
of 2/5 excites the imaginary parts associated with these fixed points, thus inducing inertial 
oscillations which are unphysical for the case of pure shear flow. 

An overview of the performance of the models can be gleaned from Table 1 which com- 
pares the predicted equilibrium values with the most recent experimental data of Tavoularis 
and Karnik 7 for homogeneous shear flow (this data constitutes a mean over the stronger 
shear rate cases). Here (•)„, denotes the equilibrium value obtained in the limit as t —* 00. 
Several observations concerning Table 1 are noteworthy: 

(a) The revised RNG K — e model yields substantially better results than the old version 
of the model and is, on balance, better than the standard K — e model. This appears to 
explain why the models performed as they did in Figure 1 relative to the LES results. 

(b) The only deficiency in the predictions of the new RNG K—t model for homogeneous shear 
flow are in the values of the normal components of the anisotropy tensor - a shortcoming of 
any model based on an isotropic eddy viscosity. However, the RNG based anisotropic eddy 
viscosity model of Rubinstein and Barton 10 - which predicts (fcu)oo = 0.260 and (622)00 = 


4 



—0.196 for the normal anisotropies in homogeneous shear flow - alleviates this deficiency to 
a large extent. 

(c) The RNG second-order closure model does perform somewhat better with the new value 
of C c \ (the model now predicts a weak exponential time growth of K* whereas the old 
version of the model predicted a power law growth, with {^SK / t)oo = oo, due to the close 
proximity of C c \ to the bifurcation point C e i = 1). However, this preliminary model still 
performs weakly in comparison to the more commonly used second-order closures such as 
the Launder, Reece, and Rodi model. The deficiency in this model is traced to the rapid 
part of the pressure-strain correlation which is 0(1) instead of 0(b) in the anisotropy tensor. 
In fact, the deviation of O 2 from 0.4 to 2/15 results from the model trying to compensate 
for the truncated 0(b) terms 11 (interestingly enough, if O 2 is set to 0.4 in Eq. (11), the 
predictions of the model deteriorate substantially). Hence, we have little doubt that if the 
RNG based calculation is extended to include the O(b) terms, the resulting model would 
perform quite well in comparison to other second-order closures. 

In conclusion, with the revised coefficients proposed by Yakhot and Smith 3 , the RNG 
K — e model now performs well in homogeneous shear flow — particularly when the RNG 
based anisotropic eddy viscosity of Rubinstein and Barton 10 is used. The RNG second-order 
closure model needs further development, however. It would appear that an extension of 
the rapid pressure-strain correlation to include terms of O(b) would resolve the remaining 
deficiency in this model. Consequently, our current assessment of RNG based turbulence 
models is now more optimistic than reported earlier. 
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(&n)oo 

(^12)00 

(^22)00 

(¥)„ 

New RNG 
K-e Model 

0 

-0.185 

0 

4.38 

Old RNG 
K-e Model 

0 

-0.489 

0 

11.70 

Standard 
K — e Model 

0 

-0.217 

0 

4.82 

New RNG 
Second- Order 
Closure 

0.489 

-0.091 

-0.244 

8.94 

Old RNG 

Second- Order 
Closure 

0.533 

0 

-0.267 

OO 

LRR Model 

0.193 

-0.185 

-0.096 

5.65 

Experimental 

Data 

0.21 

-0.16 

-0.13 

4.8 


Table 1. Comparison of the equilibrium values of the various models with the experimental 
data of Tavoularis and Karnik 7 on homogeneous shear flow. 
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